Two nonmagnetic impurities in the DSC and DDW state of the cuprate 
superconductors as a probe for the pseudogap 
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The quantum interference between two nonmagnetic impurities is studied numerically in both the 
d-wave superconducting (DSC) and the d-density wave (DDW) state. In all calculations we include 
the tunnelling through excited states from the Cu02 planes to the BiO layer probed by the STM tip. 
Compared to the single impurity case, a systematic study of the modulations of the two-impurity 
local density of states can distinguish between the DSC or DDW states. This is important if the 
origin of the pseudogap phase is caused by preformed pairs or DDW order. Furthermore, in the 
DSC state the study of the LDOS around two nonmagnetic impurities provide further tests for the 
potential scattering model versus more strongly correlated models. 
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I. INTRODUCTION 



The study of magnetic and nonmagnetic impurities in 
the Cu0 2 planes of the High-T c superconductors is far 
from settled. Experimentally, the local density of states 
(LDOS) measured by scanning tunneling microscopy 
(STM) in Bi 2 Sr 2 CaCu 2 8+5 (BSCCO) around a non- 
magnetic impurity such as Zn displays a sharp peak 
close to the Fermi level on the impurity site and a second 
maximum on the next-nearest neighbor sites 1,2 . The- 
oretically, the question remains whether a traditional 
potential scattering formalism 3,4 or more strongly corre- 
lated models 5 are needed to describe the impurity effects. 
Though still a subject of controversy, it was recently 
shown that at least for weak impurities a potential scat- 
tering scenario qualitatively agrees with the measured 
results for optimally doped BSCCO 6-10 . Furthermore, it 
was shown by Martin et ai. 11 that both the energetics and 
the spatial dependence of the resonance state around a 
strong potential scatterer (e.g. Zn) can be accounted for 
by including the tunnelling (the filter) through excited 
states from the Cu0 2 planes to the top BiO layer probed 
by the STM tip 12 . There is also evidence from nuclear 
magnetic resonance (NMR) measurements that magnetic 
moments are induced around nonmagnetic impurities 13 . 
In this paper we assume, however, that the large poten- 
tial scattering off the impurity site itself is dominating 
the final LDOS. 

Recently the experimental ability to manipulate the 
positions of surface impurities has increased the inter- 
est in quantum interference phenomena between mul- 
tiple impurities. This includes the physics of quan- 
tum mirages 14 and various multiple impurity effects in 
superconductors 15-18 . For example, it was shown in Ref. 
17 that impurity interference can be utilized as a sensitive 
probe for the gap symmetry of exotic superconductors. 
Motivated by the experimental progress we compare the 
expected LDOS around one and two strong nonmag- 
netic impurities in either the d-wave superconducting 



(DSC) or the d-density wave (DDW) state. Though 
still controversial we include the filter effect in all the 
calculations presented below. As has become clear only 
recently 15-17 , we stress that the probed impurities need 
be well separated (10-50 lattice constants) from other 
possible defects. 

The DDW state was recently proposed as a model for 
the pseudo-gap state of the cuprates 19 . Any difference 
in the impurity modified LDOS between the DSC and 
DDW states may reveal the hidden DDW order and dis- 
tinguish between the scenario of preformed pairs versus 
static staggered orbital currents as the origin for the 
pseudo-gap state 20-22 . Recently, there has been sev- 
eral other proposals to probe the DDW order in the 
cuprates 23,24 . 



II. MODEL 

In this section we briefly discuss the models for the 
DSC and DDW states and how to calculate the LDOS 
around several impurities. The BCS Greens function 
G°(k, iu> n ) for the unperturbed d-wave superconductor 
is given by 



G°(k, iuj n ) = [iu n t - £(k)f 3 - A(k)fi]" 



(1) 



where t„ denotes the Pauli matrices in Nambu space, tq 
being the 2x2 identity matrix, £(k) the quasi-particle 
dispersion, and u) n is a Matsubara frequency. For a 
system with d 2; 2_ 1) 2-wave pairing symmetry, A(k) = 
^ (cos(fc x ) - cos(fc„)). 

In the DDW state the mean-field Hamiltonian is given 
by 19 



kcr 



+ i 



kcr 



D(k)cl a Cu+Qa 



(2) 



where c^ a creates an electron with momentum k and spin 
a, Q = (7r,7r) and D(k) = (cos(k x ) —cos(k y )). Bc- 
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low, Ao = Do = 50meV and the lattice constant is set 
to unity. The large value of the gap corresponds roughly 
to the experimentally measured maximum gap in the un- 
dcrdoped regime of BSCCO. 

The Greens function for the clean DDW state is given by 



G°(k,iw n ) = 



iw n -<£(k + Q) _i£)(k) \ 
iD(k) iu^ - g(k) J 

(iw„ - £(k))(iw„ - £(k + Q)) - £)(k)2 ' 



(3) 



Performing the Fourier transform, G°(rj, rj, iu n ) = 
Ekk'G' (k,k',iw„)e lkr '- lk ' - r \ of the Greens function 
with reference to the entries of Eqn. (3) gives 

G°(r,,r,^ n ) =^[G? 1 (k,iw n )+G? 2 (k,ia; n )e- i ^ + 

k 

G° 21 (k,iuj n )e iCi ^ +G a 22 (k,iLu n y Q - (r *- r i ) }e lk - < - r *- r i\ (4) 

with the sum extending over the reduced Brillouin zone. 
The presence of scalar impurities is modelled by the fol- 
lowing delta-function potentials 



{i}a 



Uifli 



(5) 



where hi a is the density operator on site i. Here {i} de- 
notes the set of lattice sites hosting the impurities and 
Ui is the strength of the corresponding effective poten- 
tial. In this article all the presented results arise from 
impurities modelled by a potential, U = —15t, corre- 
sponding to -4.5cV. In the DSC state this U gener- 
ates resonances at a few meV for a single nonmagnetic 
impurity 1,2 ' 11 . The large scale of this potential renders 
the effects on the LDOS from corrections to other energy 
scales around the impurity site less important. For in- 
stance, we have checked that gap suppression near the 
impurity or slightly larger spatial extension of the im- 
purity does not qualitatively affect the results reported 
below. In general these effects tend to push the reso- 
nances slightly further towards zero bias. We have also 
performed calculations (not shown here) with other val- 
ues of U and comment on the results below. 
The full Greens function G(r, u) in the presence of the 
impurities can be obtained by solving the real-space 
Gorkov-Dyson equation 



G{lo) = G°(lu) (/ - H mt G\io)^j 



(6) 



The size of the matrices in this equation depends on the 
number of impurities and the dimension of the Nambu 
space. We have previously utilized this method to study 
the electronic structure around impurities 17 and vortices 
that operate as pinning centers of surrounding stripes 27 . 
This method is identical to the traditional T-matrix for- 
malism. However, for a numerical study of several impu- 
rities at arbitrary positions we find it easier to solve Eqn. 



6 directly. 

The 2D Fourier transform of the clean Greens function 
G° (k, lj) is performed numerically by dividing the first 
Brillouin zone into a 800 x 800 lattice and introducing a 
quasi-particle energy broadening of 6 = lmeV with 5 de- 
fined from the analytic continuation iu n — > lo + iS. The 
differential tunnelling conductance is proportional to the 
LDOS which is determined from the imaginary part of 
the full Greens function. 

So far nothing has been said about the form the band- 
structure. It is still controversial which quasi-particle 
energy applies to the DSC and DDW states 22 ' 25,26 . The 
expression for £(k) is important since it will influence the 
final LDOS around the impurities. We illustrate this in 
the following by studying two generic band structures: 
the nested situation, and a t-t' band believed to be rel- 
evant for BSCCO around 10% hole doping. With the 
notation £(k) = e(k) — fi, and 

e (k) = — 2t (cos(k x ) + cos ( k y ) ) — At ' cos ( k x ) cos ( k y ) , (7) 

t (t') refers to the nearest (next-nearest) neighbor hop- 
ping integral and /i is the chemical potential. The nested 
situation corresponds to t' = /j, = 0.0 while the param- 
eters for the 10% hole doped band are: t = 300meV, 
t' = —0M and /j, = — 0.9i. These parameters corre- 
spond to the ones previously studied for a single impurity 
by Morr 22 . As discussed in Ref. 22 there are physi- 
cal reasons to expect the nested band to be relevant for 
the DDW state and the t-t' band for the DSC phase. 
However, recent photoemission measurements on LSCO 
by Zhou et al. 26 observed a Fermi surface consisting of 
straight lines connecting the antinodal regions which may 
indicate that the nested band is more relevant for impu- 
rity studies in LSCO. Thus we find it important for study 
both cases below. 

In the results presented we include the LDOS 
filter 11 . This effect modifies the LDOS, p(r,u) = 
X) n IV'ra( r )| 2 ^( w ~ e n), by including the four nearest Cu 
neighbors in the underlying Cu02 layer, ip n { r ) — ► 
■0„(r + e x ) + tjj n (r - e x ) - ip n (r + e y ) - -0„(r - e y ). 
Here denote the unit vectors in the Cu02 plane. It 
is important to keep in mind that the filtering effect is 
still controversial. However, determining experimentally 
the interference effects around two impurities in the DSC 
state may help resolve the relevance of the filter. 



III. RESULTS 

A. single impurity 

Before studying the two impurity interference effects it 
is worthwhile to briefly review the single impurity LDOS 
in the DSC and DDW states and discuss the influence of 
the filter. Without the tunneling filter we find full agree- 
ment with previously published results 7, 15, 20-22 . We will 



2 



see that a single impurity is not a good probe for distin- 
guishing between these two states. 

In the DDW phase one can utilize Eqn. (4) and (6) 
to calculate the full Greens function G(rj, r^, iw n ) = 
G°(r l -r J ,iuj n )+G a (r il iu n )T(iu) n )G°(-r jl iLo n ) with the 
T-matrix given by, T(iu n ) = U[l - UG°(0, iw n )] _1 . The 
single resonance condition, 1 = UKc[G°(0,co)], has been 
previously studied for the DDW state without the filter- 
ing effect 20-22 . It is well known that the resulting LDOS 
strongly depends on the band structure. In Fig. la we 
plot the DOS in the clean DDW state for the nested 
and the t-t' band without the filter. Even though the 
above resonance condition is satisfied at certain energies 
for the t-t' band, we expect the large value of the DOS at 
all frequencies to overdamp the impurity peaks. This is 
contrary to the nested situation where a sharp impurity 
resonance is allowed to appear in the gap. This is clearly 
verified in Fig. lb(c) which depicts the LDOS for the 
nested(i-i') set of band parameters including the filter. 
The peaks in Fig. lc are not impurity resonances (note 
scale), which are overdamped, but simply the shifted 
DDW gap edges. The impurity can only slightly mod- 
ify the amplitude of these gap edges. We note that it is 
t' which causes the impurity resonances to be strongly 
overdamped. When t' = 0, fi ^ the density of states 
always vanishes at minus the chemical potential 21 allow- 
ing a well-defined resonance peak to appear. 
As is evident from Fig. lb the most important influ- 
ence of the filter is to shift the LDOS maximum from the 
nearest neighbors to the impurity site and induce a sec- 
ond maximum on the next-nearest neighbor sites 28 . This 
weight redistribution is identical to the situation in the 
superconducting phase 11 . 

In the DSC state, the clean DOS is plotted in Fig. 2a 
for both the nested and the t-t' band. By comparison 
to Fig. la we see the well known result that the nested 
DOS is identical for the clean DDW and DSC phase. In- 
deed this motivated the original studies of single impurity 
resonances in the DSC versus DDW states 20-22 . The 
single impurity resonance condition in the DSC phase, 
1 = f7Re[G°(0, ±o»)], generates peaks at positive and 
negative energies around a single nonmagnetic impurity. 
However, the majority of the quasi-particle weight may 
reside on only one of these resonances 17 . 



(b) 


(c) 
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FIG. 1. DOS (arb. units) as a function of energy (meV) in 
the DDW state: (a) for the clean system with nested (solid) or 
t-t' (dashed) band, (b) DOS at (0,0) (solid), (1,0) (dashed), 
and (1, 1) (dash-dotted) for a nested band with the impurity 
at (0,0). (c) same as (b) but for a t-t' band. 
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FIG. 2. Same as Fig. 1 but for the DSC state. 



It is evident from both Fig. 2b and Fig. 2c that indeed 
only one resonance has weight. This is contrary to the 
situation without the filter 22 . Thus by comparing Fig. 
lb to Fig. 2b (or 2c) the result is two almost identical 
figures. Therefore, since no qualitative difference is guar- 
anteed to exist the single nonmagentic impurity cannot 
easily distinguish the DSC and DDW phases . However, 
as shown below, the interference between several impuri- 
ties can be utilized to tune the amplitude of the potential 
resonances and thus clearly distinguish the phases. 
The impurity LDOS plotted in Fig. 1 and Fig. 2 was 
for U = —15t. Though of less experimental relevance, 
we briefly mention another difference between the DSC 
and DDW states. This relates to the fate of the reso- 
nance in the unitary limit, U — > oo: for the DDW phase 
the resonance energy approaches minus the chemical po- 
tential, u> — —fi, whereas it approaches the Fermi level 
in the DSC phase (except for a small residual energy 
shift caused by a possible particle-hole asymmetry 29 ). 
The different resonance energy (as U — ► oo) arises from 
the way the chemical pot ential enters th e bands of the 
clean DDW {E±(k) = | ^(k) 2 + L>(k) 2 ± fj,\) and DSC 
(E ± (k) = v /(|e(k)|± A i) 2 + A(k) 2 ) states 20 ' 21 . 



B. two impurities, nested band 



In general when several impurities are in close proxim- 
ity the resonances split, and one expects to see additional 
peaks in the density of states. The evolution of the LDOS 
as a function of distance and angular orientation between 
two nonmagnetic impurities in the DSC state has been 
already studied by several authors 15 17 . In the follow- 
ing we elaborate on this work by a numerical study of 
the LDOS including the filtering effect and study for the 
first time the quantum interference between two strong 
nonmagnetic impurities in the DDW state. In the super- 
conducting phase Fig. 3a shows the resulting LDOS for 
the nested band when one impurity is fixed at the origin 
(0, 0) while the other is moved out along a crystal axis to 
(10, 0). In Fig. 3b the impurities are fixed at (—1, 0) and 
(+1,0) while the STM tip is moved from (0,0) to (8,0). 
As seen from both figures there are strong variations in 
the LDOS in agreement with previous studies without 
the extra tunnelling effect 15 ' 17 . 
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FIG. 3. (a) DOS at (0, 0) as a function of the distance 
between two nonmagnetic impurities (here: t' = p, = 0.0). 
One impurity is fixed at (0, 0) while the other moves from 
(0,0) (top) to (10,0) (bottom), (b) the impurities are fixed 
at (±1,0) while the STM tip is moved from (0,0) (top) to 
(8,0) (bottom). The difference between each scan is 0.2 lat- 
tice constants and the graphs are off-set for clarity. 
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FIG. 4. Same as Fig. 3 but for the d-density wave state. 

The number of apparent resonances, their energetic posi- 
tion and width strongly depend on the impurity configu- 



ration and the position of the STM tip. In particular, for 
certain impurity separations the resonances completely 
disappear. In Fig. 4 we show the LDOS for the same 
impurity and STM positions as in Fig. 3 but for the 
DDW state. Clearly, strong quantum interference be- 
tween the two nonmagnetic impurities also exists in this 
state. However, by comparison with Fig 3 it is evident 
that the additional resonance states in the DSC allows 
one to distinguish this from the DDW phase. We have 
performed identical calculations to the ones presented in 
Fig. 3-4 for other (but still large) values of the scatter- 
ing potential U, and always find qualitatively the same 
interference pattern. 

As mentioned above, the resonances split when two im- 
purities arc in close proximity. It is therefore nontriv- 
ial that only a single, nondispersive peak is seen in e.g. 
Fig. 4b. This is closely connected to the particular 
STM scan and one may worry about the robustness of 
this result. However, we always find that whenever the 
impurity positions are invariant under mirror reflection 
through the STM scan line, only a single nondispersive 
peak remains 30 in the DDW state. Importantly, for these 
same configurations we find the alternating double peak 
structure (similar to Fig. 3b, 5b) to be a robust feature in 
the superconducting phase. Furthermore, as expected for 
a d-wave gap 17 , we find (not shown) that the quantum 
interference patterns are longer ranged along the nodal 
directions than along the Cu-0 bonds. 
As expected from the discussion of the single impurity in 
the DDW state, we end this section by noting that when 
if = 0.0, fi 0.0 the interference pattern is identical to 
that shown in Fig. 4 except for a shifted (by —fx) energy 
range. 



C. two impurities, t-t' band 

We now turn to the quasi-particle dispersion given by 
Eqn. (7) with t' = -0.3t, Li = -0.9t. In this case we 
know from Fig. la and Fig. 2a that the clean DOS 
is clearly different in the DDW and DSC states. This 
section serves as an illustration of the importance of 
the quasi-particle dispersion in the final LDOS. Fig. 5 
shows the LDOS in the superconducting phase from the 
same STM and impurity positions as Fig. 3. It is clear 
that again the strong interference between the impurity 
wavefunctions survive the filtering effect and pose new 
constraints on the potential scattering scenario versus 
more strongly correlated models 5 . We note that despite 
the very different band structure used to calculate the 
LDOS in Fig. 3 and Fig. 5, the overall evolution of the 
resonances is quite similar except that the apparent res- 
onances are shifted to higher energies for the t-t' band. 
As mentioned above, it has been previously suggested 
that the nested (t-t 1 ) band is appropriate for the DDW 
(DSC) state 22 . In that case we need compare Fig. 4 and 
Fig. 5. As opposed to the single impurity LDOS, the 
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configuration in Fig. 4b and Fig. 5b again allows one to 
distinguish the DDW and DSC states by the number of 
resonance peaks. This is contrary to Fig. 4a and Fig. 5a 
which are remarkably similar. 
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FIG. 5. Same as Fig. 3 but for t' = -0.3t and \i = -0.9t. 
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FIG. 6. Same as Fig. 5 but for the d-density wave state. 
Note the energy range. The peaks evident in these scans are 
the shifted DDW gap edges (not impurity resonances) . 



In the DDW phase we know from the single impurity 
case that the current choice of band parameters leads 
to strongly overdamped impurity resonances (Fig. lc). 
However, for completeness we show the calculated STM 
scans in Fig. 6. As expected the quantum interference 
is weak and causes only minor changes in the DDW gap 
edges. Furthermore, the LDOS shown in Fig. 6 changes 
only slightly upon varying U or the impurity positions. 



IV. CONCLUSION 

In summary we have shown that a systematic STM 
study around two nonmagnetic impurities can clearly 
distinguish the DSC and DDW phases. In particular, 
we suggest to perform STM scans with the positions of 
the impurities being invariant under a mirror reflection 
through the scan line. Even for the nested band, where 
the clean and the single impurity LDOS are not a good 
probe for the underlying state, this situation provides a 
robust test for DSC versus DDW order. 
The impurities are modelled as potential scatterers and 
the results pose further tests on this approach. An impor- 
tant question remains whether phase fluctuations present 
about T c in the pseudo-gap state are strong enough to 
wash out the interference patterns. This will be discussed 
in a future publication 31 . It would also be interesting 
to study similar multiple impurity interference effects 
within other pseudo-gap models and within other pro- 
posed scenarios for the resonances around nonmagnetic 
impurities in d-wave superconductors. In particular, 
within models explaining the single impurity LDOS as a 
Kondo resonance arising from a confined spinon 5,32 , one 
may expect more novel changes as the distance between 
two nonmagnetic impurities is decreased. This is because 
the cost of frustrated dimers decrease in this limit mak- 
ing it unfavorable to break another dimer, and hence no 
spin is expected near the nonmagnetic impurities. 
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Hedegard, J. Paaske, S. Sachdev, and J.-X. Zhu. This 
work is supported by the Danish Technical Research 
Council via the Framework Programme on Superconduc- 
tivity. 
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